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We calculate the strength of the electroweak phase transition in a supersymmetric model with 
four chiral generations. The additional chiral fermions (and scalar partners) lower the critical 
temperature and thus strengthen the first-order phase transition. The scalar partners stabilize the 
potential, leading to an effective theory that is bounded from below. We identify the ensemble 
of parameters where (j> c /T c ^ 1 simultaneous with obtaining a large enough Higgs mass. Our 
calculations focus on a subset of the full four generational supersymmetric parameter space: We 
take the pseudoscalar heavy, tan/3 = 1, and neglect all subleading contributions to the effective 
potential. We find that the region of parameter space with a strong first-order phase transition 
requires m^i/m q i J$ 1.1 while the constraint on the lightest Higgs mass requires m q i/m q i > 1 with 
m q i > 300 GeV. We are led to an intriguing prediction of quarks and squarks just beyond the 
current Tevatron direct search limits that are poised to be discovered quickly at the LHC. 



I. INTRODUCTION 

The origin of the matter asymmetry is a deep mystery 
that remains unsolved. Conditions that can lead to a dy- 
namical asymmetry between baryons and anti-baryons 
were articulated years ago by Sakharov pQ: baryon num- 
ber violation, C and CP violation, and a departure from 
thermal equilibrium. All three conditions are satisfied by 
the Standard Model as it passes through the electroweak 
phase transition. But, the CP violation is too small [2], 
and the phase transition is not strongly first-order (e.g., 
[3J HJ |5J |5J |7]), given the direct search bounds on the 
Higgs boson from LEP. New physics with large CP vi- 
olation is trivial to introduce into the model; weak scale 
supersymmetry is an obvious example (care needs only 
to be taken to ensure that induced electric dipolc mo- 
ments are within the experimental bounds) . Even with a 
new source of CP violation, if the phase transition is not 
strong enough first-order, any generated baryon asymme- 
try will be washed out [5]. New physics that enhances the 
first-order phase transition, however, is generally much 
more tricky to achieve. 

In the early 1990s it was realized that the electroweak 
phase transition could be enhanced by modifying the 
cubic coupling of the finite-temperature effective poten- 
tial [3]- Nontrivial modifications of the cubic coupling 
could arise from additional scalars with order one cou- 
plings to the Higgs. In the minimal supersymmetric stan- 
dard model (MSSM), the scalar superpartners to the top 
quarks (stops) can play precisely this role [5J QIJJ E] • It 
has long been advocated that the region of MSSM pa- 
rameter space with a light stop (and a light Higgs) can 
yield a strong enough phase transition. Unfortunately, 
the combination of direct searches for the lightest Higgs 
boson and direct searches for stops have virtually ruled 
out this possibility. The remaining parameter space |12j 
requires a large hierarchy between the left-handed and 
right-handed stops to ensure the Higgs satisfies the LEP 
bound. 

Methods to strengthen the first-order phase transition 
beyond the SM and MSSM are now widely discussed [T3J 



HH HS1 H6l HZl [HI [HI [20I [21]. Several of these ideas 
add a singlet field, such as in the NMSSM or nMSSM. 
Another related idea is to simply cutoff the SM at a low 
scale, adding the effects of higher dimensional operators 
[22] (which can be UV completed by integrating out a 
singlet) . 

Yet another interesting possibility, and the one we will 
focus on in this paper, is to add more particles with mod- 
estly strong couplings to the Higgs. This was proposed in 
[T7] ; the additional particles have quantum numbers such 
that they mix with the MSSM charginos and neutralinos. 
Heavy particles that receive their mass entirely or dom- 
inantly from electroweak symmetry breaking can have a 
substantial impact on the electroweak phase transition. 
In this paper we consider a modification to the MSSM 
similar in spirit to [17] to enhance the phase transition. 
Namely, we add a fourth generation of particles (and 
sparticles) to the MSSM. Larger couplings to the Higgs 
are automatic simply due to the direct search bounds 
from LEP and Tevatron on fourth generation fermions. 

That the electroweak phase transition could be en- 
hanced in a four generation supersymmetric model was 
considered before in [23] , They performed an interest- 
ing numerical study that also found that the electroweak 
phase transition can be enhanced when rriqi jm q i is not 
much larger than 1. However, the present limits on fourth 
generation quark masses rule out their parameter space, 
and moreover, they allowed tan/3 to far exceed 1, imply- 
ing that the b' Yukawa coupling was nonperturbative. In 
our analysis, we first systematically analyze the origin of 
the contributions that allow the phase transition to be- 
come first-order. This allows us to make a clear distinc- 
tion between supersymmetric and non-supersymmetric 
theories with heavy chiral fermions. We then identify 
the viable region of parameter space where the Yukawa 
couplings are under control (tan/3 ~ 1) and all other 
bounds are satisfied. 

A fourth generation has historically been thought be 
strongly disfavored by the absence of flavor mixing, the 
Z — > vv constraint, and electroweak precision data. All 
of these objections can be straightforwardly overcome, 
as was recently emphasized in [24] . Below we summarize 



these results in the form of the parameter space that is 
allowed. It is interesting that the constraints from elec- 
troweak precision data can be overcome without or with 
an electroweak scale Majorana mass for the fourth gen- 
eration right-handed neutrinos. If a Majorana mass does 
indeed exist , lepton number is violated at the electroweak 
scale, and thus scenarios of baryogenesis that rely on an 
earlier generation of B — L number (such as leptogene- 
sis) do not work here [251 US]- Electroweak baryogenesis 
is one of the few mechanisms not sensitive to this source 
of lepton number violation, and thus becomes even more 
interesting to study. 
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II. SETUP 

We consider a low energy supersymmetric theory with 
a fourth chiral generation of matter, the "4MSSM" (for 
an early discussion, see |27J. A fourth chiral generation 
of matter does affect electroweak precision observables. 
One of the main results of [24] was to show that it can 
be made completely consistent with electroweak precision 
data so long as there are modest mass splittings between 
the isospin partners in the quark and lepton doublets. 
This splitting causes a modest reduction in the positive 
contribution to S, simultaneous with a positive contri- 
bution to T, allowing ensembles of parameters that are 
within the 68% CL ellipse of the LEP Electroweak Work- 
ing Group. For example, the mass spectrum m u i = 100 
GeV, m t > = 155 GeV, m t > = 310 GeV, m v = 260 GeV, 
and rrih = 115 GeV is perfectly acceptable. There is 
strong sensitivity to the mass differences with only mild 
sensitivity to the overall scale of the particles. We will 
present results for both the electroweak preferred ratio 
frit' /rn,b' ~ 1-2 as well as mf = my for comparison. 

With supersymmetry, there are additional contribu- 
tions to electroweak precision observables from super- 
partners (e.g., [28). We have not included these con- 
tributions to optimize the parameter set to match elec- 
troweak data, simply because many more parameters en- 
ter the fit that can be freely adjusted without affect- 
ing our results for the electroweak phase transition. We 
therefore take fourth generation Yukawa couplings con- 
sistent with [53] and take the scalar partner masses to be 
degenerate, eliminating this potential additional contri- 
bution to isospin violation. 

We neglect all sub-leading contributions to the zero- 
temperature and finite temperature effective potential. 
Sub-leading here refers to couplings smaller than about 
1. We retain the contributions from gauge bosons. But 
we neglect light fermions (u,d,c,s,b,e,/x,r), Higgs bosons 
(the quartic is small), and all superpartners other than 
i'1,2, &'i,2, *i,2- We also neglect contributions from fourth 
generation leptons because the number of degrees of free- 
dom per particle is only 1/3 that of quarks and the 
bounds on the mass from the non-observation in experi- 
ment are much weaker than for quarks. 



FIG. 1: An estimate of the cutoff scale of the 4MSSM as a 
function of the fourth generation fermion mass by running 
the one-loop RGEs of the quark Yukawas up to where they 
encounter a Landau pole (j/y/ ~ 47r). 

III. SUPERSYMMETRY WITH tan/3 = 1 

In the limit tan (3 — > 1, several aspects of supersymme- 
try drastically simplify. From the definition of tan (3 = 
(H u ) I '(Hj) , we see the vevs are equal, v u — Vd = v/y2, 
where v = 246 GeV, and the Yukawa couplings 



are the same for the up-type and down-type fermions. 
The fourth generation quarks t' and b' have large Yukawa 
couplings, 
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where 260 GeV is shown for convenience in comparison 
to the (approximate) present direct search bounds from 
the Tevatron [20J 130] ■ Note that these Yukawa couplings 
are a factor of y2 larger than in a non-supersymmetric 
model, since the t' and b' acquire their mass only through 
couplings to the up-type and down-type Higgs, respec- 
tively. 

If tan/3 7^ 1, either the up-type or down- type Higgs 
vev is reduced, and thus to hold the masses of the 
fermions fixed, either y t i or y^ must increase. Setting 
tan (3 = 1 allows the largest possible physical fourth gen- 
eration fermion masses with the smallest Yukawa cou- 
plings. Since y/ cannot be arbitrarily large for pertur- 
bative calculations at the weak scale to be valid, the pa- 
rameter choice tan (3—1 really just maximizes the cut- 
off scale of the model. Even with this adjustment, the 
cutoff scale is low. This can be estimated by running the 
one-loop renormalization group equations for the Yukawa 
couplings up to ~ 4-7T. We show the scale of the Landau 
pole as function of fermion mass in Fig. [I] Note that 
requiring y'j/^Air) < 1, implies yf < 3.5, corresponding 
to mj < 450 GeV; we will not consider fermion masses 
that much exceed this value. 
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In the limit tan (3 — > 1 the Higgs sector also drastically 
simplifies. The tree-level potential in the MSSM with 
tan/3 = 1 is [35 
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Expanding the neutral components as 
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where the a rotation matrix diagonalizes the Higgs fluc- 
tuations (h, H) into mass eigenstates. In the limit 
that the second Higgs doublet "decouples" (m A o^ H _ H ± » 
m/j), the mixing angle a — > —(3, and thus the lightest 
Higgs is simply h = (iJ^ + H® — \/2v). In this limit the 
tree-level Higgs potential vanishes, since h corresponds 
to the excitation of a Z3-flat direction. 

Since electroweak precision data prefers m t i jmy ~ 1.2, 
this could be arranged either by adjusting just these 
two Yukawa couplings yv jyv = 1.2 or instead adjusting 
tan/3 = 1.2. These two scenarios are nearly equivalent 
for our purposes, and so we choose to set tan (3=1. The 
alternative, tan/3 = 1.2, would give a tree-level contri- 
bution to the Higgs potential. The contribution to the 
(mass) 2 is however just 0.03A/§. As we will see, the one- 
loop radiative corrections from quarks and squarks will 
be far larger than this, so it is safe to completely neglect 
tree- level contributions even if tan /3 were allowed to vary 
slightly from 1. 

In addition to taking tan/3 = 1, we also choose super- 
symmetric parameters such that the mass eigenstates of 
#1,2 and &'i,2 correspond to the gauge eigenstates t'z,R 
and b 1 l,r- This is done purely to simplify our calcula- 
tion. It is a rather conservative approximation, since it 
is well known that increasing the off-diagonal contribu- 
tion to the squark mass matrix leads to an enhancement 
in the one-loop contribution to the Higgs mass (e.g., see 
[5"2"]). We expect that the parameter space with a strong 
first-order phase transition will enlarge as this restric- 
tion is relaxed. Note that since the off-diagonal left-right 
contribution to the up-type and down-type squark mass 
matrix is equal to nif(Af — fJt) (where again, tan/3 =1), 
this simplification corresponds to the specific parameter 
choice Af = /i. 

Finally, as we discussed above, the Higgs potential sim- 
plifies in the limit m A o,m H ±,mjjo 3> m^. This is a 
common assumption in the electroweak phase transition 
literature: The operational advantage is that the low en- 
ergy theory is effectively a one-Higgs-doublet model that 
is drastically simpler to analyze at finite temperature. 



IV. ONE-LOOP EFFECTIVE POTENTIAL 

In the 4MSSM with tan (3=1, loop corrections entirely 
determine the Higgs potential. We are interested in the 



loop corrections to just the scalar fluctuation cj) = (h + 
v) At one-loop the effective potential for the Higgs 
is determined from the Coleman- Weinberg potential 
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where M^ifi) are the field-dependent masses and fj, is the 
renormalization scale (MS scheme). We generally use 
Mi to refer to ^-dependent (and temperature-dependent) 
masses and mi to refer to Mi(v) at zero temperature. The 
c/s are constants corresponding to 5/6 for gauge bosons 
and 3/2 for fermions and scalars. The degeneracies per 
particle are n q = —12 (for each q — t, t' , b'), riq L = n§ R = 
6, n WT = 4, u Zt = 2, u Wl = 2, n Zh = 1. 

Expanding the effective potential as given above, ev- 
idently the minimum is not necessarily located at the 
proper electroweak breaking scale v = 246 GeV. This 
is easily remedied by imposing a renormalization con- 
dition on the mass parameter such that the minimum 
is enforced to be at v. This amounts to adding the in- 
dependent contribution to the effective potential, 
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The masses used in the effective potential are MS 
masses that differ from the physical (pole) masses 
through finite and log-dependent corrections. The run- 
ning fcrmion masses are given at one-loop by 
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Since the fourth generation fermions overwhelmingly 
dominate the contributions to the (zero- and finite- 
temperature) effective potential, we take /i = ^/rnt^ffiy, 
i.e., the scale of the largest electroweak breaking masses 
in the problem. This tends to minimize the higher order 
corrections to the potential, though are calculations are 
not particularly sensitive to the precise choice of renor- 
malization scale. 

The running scalar masses also differ from their phys- 
ical pole masses through one-loop corrections depending 
on not only the gluon but also gluino diagrams |28j . This 
correction is generally numerically smaller than the cor- 
rection to the fermion mass, typically less than a few %. 
Moreover, since the correction is gluino mass-dependent, 
relating the pole mass to the running mass requires spec- 
ifying an otherwise unfixed parameter in our model. We 
choose instead to simply take wj polc = mj(mj), thus 
neglecting the difference between the pole and running 
mass for the squarks. 
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FIG. 2: An illustration of the effect of adding one additional heavy fermion that obtains a mass of 400 GeV entirely from 
electroweak symmetry breaking. The figure (a) on the left shows the effective potential at the critical temperature T c (solid 
line) in the 4MSSM with m t i = my = 400 GeV and m,-/ /m q i = 1.05. Adding an additional heavy fermion (with n = — 1) to the 
effective potential is shown (dashed- line), except that only the finite-temperature contribution, Vri, is included in (a). Figure 
(b) on the right shows the the effect of including just the finite-temperature contribution (dashed-line) , identical to Figure (a), 
and then the effect of including both the finite-temperature as well as the zero-temperature contribution Vi (solid-line). The 
net effect shown in Figure (b) solid-line is that the global minimum at <f> ~ 4> c decreases and thus T c increases. 



V. FINITE TEMPERATURE ONE-LOOP 
EFFECTIVE POTENTIAL 

The finite-temperature contributions to the effective 
potential are [3 ES G3 EE E] 



evaluate the thermal functions Jb,f numerically. The 
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The field-dependent fermion and scalar masses are 
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Explicit expressions for the thermal masses of the SM 
gauge bosons can be found in, e.g., Ref. [33]. The finite- 
temperature contributions depend on the thermal func- 
tions 
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Often a high temperature approximation is employed to 
estimate these integrals. In our case, due to the large 
Yukawa couplings, this approximation is generally not 
appropriate. Consequently, all computations given below 



ring contribution (10 1 is only relevant for the longitudi- 



nal components of the W and Z. No contributions from 
scalars are included here since the squarks receive a con- 
tribution from soft supersymmetry breaking, and thus 
they remain heavy in the <fi —> limit. 

We have self-consistently included contributions in 
the finite-temperature effective potential the same as 
those included in the zero-temperature effective poten- 
tial. Namely, we include contributions from t',b',t, their 
superpartners t' L R , b' L R , t^ii, as well as the transverse 
and longitudinal components of W, Z. 

In the standard model, the electroweak phase tran- 
sition becomes second order when the Higgs mass ap- 
proaches 70 GeV [34] . Qualitatively, this is because 
the transverse modes of W and Z, which would drive 
a first-order phase transition in the standard model with 
a lighter Higgs mass, develop a thermal mass from non- 
perturbative effects. If the transverse thermal masses are 
large, they effectively remove the cubic term from the 
finite-temperature potential when the effective potential 
is reset to zero at (f> = 0. In our model, the first-order 
phase transition is mostly driven by squarks. In fact, 
when the W and Z are neglected in our model, <j) c and 
T c does not change significantly and the phase transition 
remains first-order. Therefore, we do not expect non- 
perturbative effects encountered in the standard model at 
larger Higgs masses to significantly affect our calculations 
of the strength of the phase transition in the 4MSSM. 



VI. EFFECTS OF NEW HEAVY PARTICLES 

The effects of heavy particles (that receive their mass 
dominantly from electroweak symmetry breaking) on the 
electroweak phase transition can be broadly character- 
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ized as follows. Consider the effective potential at T c , 
where there are two degenerate minima V e //(0,T c ) = 
V e ff((/> C ,T C ) located at <j> — and <f> = <j> c . Now add to 
this a new particle that satisfies M((f> c )/T c ^S> 1. The 
phase transition strength is modified in two ways from 
the contributions of the new particle. One is through cor- 
rections to the finite-temperature contribution; the other 
is through the zero-temperature Coleman- Weinberg po- 
tential. 



A. Finite-temperature effects 

The contributions from bosons and fermions with 
masses larger than the critical temperature, m S> T C1 can 
be characterized by how they contribute at large field val- 
ues 4> 3> T and small field values <fi <C T. At large field 
values, we can take a low temperature approximation to 
the finite-temperature effective potential. In this limit, 
the contribution from fermions or bosons becomes 



3/2 



T 4 exp 
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(15) 



where n counts the number of degrees of freedom per bo- 
son or fermion with field-dependent mass M(<fi). Clearly, 
when M((j>) ^ T, which is equivalent to <fi T (with 
order one or larger Yukawa couplings), the contribution 
to the effective potential from fermions or bosons is ex- 
ponentially suppressed. 

At small field values, we can take a high-temperature 
approximation to the thermal contribution to the effec- 
tive potential. The leading order contribution is the field- 
independent constant 
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where cg,F 
tion. 



(1,7/8) for a boson or fermion contribu- 



Thc combination of ( 15 ) and ( 16 1 imply that the intro- 



duction of a heavy fermion or boson causes a substantial 
negative shift in the potential at <p — while causing 
a negligible shift in the potential at (j> = (j> c . As an il- 
lustration, we show in Fig. [2ja) the effect of adding one 
additional heavy fermionic degree of freedom (n/ = — 1, 
for illustration) that obtains a mass of 400 GeV entirely 
from electroweak symmetry breaking. Readjusting the 
minimum V c k(4> — 0) = shifts the potential up for all 
field values, thereby removing the second minimum at 
cj) = <p c , and thus restoring electroweak symmetry. We 
must lower the temperature further in order to have the 
second minimum reappear in the effective potential with 
the new heavy fermion or boson. 



B. Zero-temperature Effects 

The second effect of heavy bosons and fermions is that 
they also modify the zero-temperature effective poten- 
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FIG. 3: This plot shows that the electroweak phase transition 
order parameter decreases if an additional heavy chiral quark 
is added to the theory without its corresponding superpartner. 



tial. Here, however, the effect of fermions and bosons 
is different. There are two contributions whose origin is 
ultimately the Coleman- Weinberg potential. One contri- 
bution is to the quartic coupling (pi, while the second 
contribution is the quadratic term M|. For smaller field 
values, i.e., <p < /i, the dominant contribution is from 
the quadratic term. Since we choose /i ~ rn q i, the log 
term drops out of Q6J) , giving an overall negative (posi- 
tive) contribution to the effective potential from fermions 
(bosons). 

The negative contribution from fermions at modest 
field values actually overpowers the effect from the finite- 
temperature contributions discussed above. This is illus- 
trated in Fig. [2jb) . The net result is that introducing a 
new heavy chiral fermion causes a decrease in 4> C /T C as 
the mass of the fermion is increased. The resulting de- 
crease in the strength of the electroweak phase transition 
with one additional chiral fermion is shown in Fig. [3] 

Adding bosonic contributions cancels the contribution 
from fermions in the Coleman- Weinberg potential. This 
cancellation is one-loop exact in the limit mj, — mf, 
i.e., no SUSY breaking contribution to the scalar mass. 
This makes it clear that we need both heavy fermions 
and scalars with equal numbers of degrees of freedom 
and similar masses to utilize the mechanism of Ref. |17j 
to lower <f> c /T c . 

Ref. [17] estimated that of order ten or more degrees of 
freedom is needed to enhance the phase transition suffi- 
ciently to achieve 4> C /T C > 1. A fourth generation quarks 
corresponds to adding 24 degrees of freedom. (We could 
have equivalently added degrees of freedom in other ways, 
such as several pairs of vector-like supersymmetric lepton 
doublets that only get mass through the Higgs mecha- 
nism. This is another interesting possibility that we will 
not explore here [35].) We have calculated the strength 
of the phase transition for a range of quark and squark 
masses. The results are shown in Fig. [4] 
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FIG. 4: Contour plot with the ratio m,-/ /m q / fixed, from top 
to bottom, as 1.05,1.1,1.15,1.2. The masses of the fourth 
generation quarks are taken to be equal. 
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FIG. 5: Contour plot of the lightest Higgs mass as a function 
of the fourth generation quark mass. Each contour corre- 
sponds to a fixed ratio m,-/ /m q i , from top to bottom, 1.2, 1.1, 
1.0, 0.95 with m t i — m b , and = ma. 



VII. LIGHTEST HIGGS MASS IN THE 4MSSM 



Given the parameter choice tan/3 = 1, the tree-level 
Higgs potential vanishes, and thus the lightest Higgs 
mass also vanishes at leading order. It is well known 
that loop corrections from splitting the masses of the top 
quark from the stops in the Coleman- Weinberg poten- 
tial provides large corrections to the tree-level value. In 
the 4MSSM, we can split not only the top and stops, 
but also split the fourth generation quarks from squarks. 
Since the one-loop contribution to the Higgs quartic cou- 
pling is proportional to yi, even a small splitting between 

/ and / has a very important effect on the Higgs mass. 
A one-loop estimate of the lightest Higgs mass in the 
4MSSM can be obtained by taking d 2 (V + Vi)/d<j} 2 at 
(j) = v. This gives our rough estimate for the Higgs mass 
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where again v — 246 GeV. In Fig. [5] we show Higgs mass 
plotted against different nif — my masses, where all 
squark masses were taken to be degenerate mjv — my — 
TOj. Each contour has the fourth generation squark-to- 
quark mass ratio, mf,/mf>, fixed to the values shown. 
Clearly, when the splitting between the fourth genera- 
tion squark and quark (MS) masses vanishes, there is an 
insufficient one-loop contribution from top/stop loops to 
raise the Higgs mass much above about 60 GeV. Never- 
theless, for even a small splitting between fourth genera- 
tion squarks and quarks, one can easily obtain a one-loop 
contribution to the Higgs mass that far exceeds the LEP 
bound so long as to/ > 300 GeV. 



VIII. RESULTS 

Combining our calculation of the phase transition with 
our calculation of the Higgs mass, we can find the al- 
lowed region in parameter space where the first-order 
phase transition is strong <fi c /T c > 1 while the Higgs mass 
satisfies the LEP bound > 115 GeV. We have com- 
puted this for the mass ratio m t '/my = 1 in Fig. [6] and 
m t > j my — 1.2 in Fig. [7] In Fig. [6]Ja) we show contours 
of increasing <p c /T c , illustrating that it is straightforward 
to obtain the strength of the phase transition to signif- 
icantly exceed (f> c /T c = 1. In Fig. |6jb) we show con- 
tours of increasing Higgs mass, illustrating that it is also 
straightforward to to obtain a lightest Higgs mass that 
significantly exceeds mu = 115 GeV. 

Note that our plots are with respect to the pole masses 
of quarks and squarks (as well as the ratio mqi/m q i). 
The quantities that enter the effective potential are MS- 
renormalized masses, which differ (as we discussed above) 
for fermion masses. Since the fermion pole mass is larger 
than its MS-renormalized counterpart by about 5%, the 
ratio of pole masses can be as small as 0.95 while the ratio 
of MS masses is still larger than one. This is why the 



fourth generation contributions to the Higgs mass (17 1 



remains positive even when the pole mass ratio rriqi / m q > 
is smaller than one. 

These results suggest that even though only one-loop 
approximations for the effective potential and the Higgs 
mass calculation were employed, we are not near any 
critical boundary, and so a more refined calculation is 
expected to only modestly adjust the parameter regions 
we have shown. For instance, there are several effects 
that can increase the Higgs mass to values higher than 
our calculation. One is moving away from the parameter 
choice Af — fi, where the off-diagonal contribution from 
squarks provides a modest increase. Another is separat- 
ing the third generation squark masses from the fourth 
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FIG. 6: The region in the quark/squark mass plane where the electroweak phase transition is strongly first-order is shown. 
The regions shown in the left- and right-hand side figures (a) and (b) are identical: the upper boundary (the solid line) is 
determined by 4>c/T c = 1 while the lower boundary is determined by m/, = 115 GeV. The dotted and dashed contours on the 
left-hand side figure (a) corresponding to <f> c /T c = 1.5,2.0 respectively. The dotted, dot-dashed and dotted contours on the 
right-hand side figure (b) correspond to the Higgs masses n%h = 150, 200, 250 GeV respectively. 



generation squark masses. The latter effect is also mod- 
est: in Fig. [5] the contour 0.95 corresponds to just the 
contribution from the third generation, and one can see 
that there is a small increase as the squark mass (equal 
to the fourth generation quark mass) increases. Simi- 
lar statements also hold for the finite-temperature effec- 
tive potential, since again our results show that there are 
model parameters where the 4MSSM model has an elec- 
troweak phase transition with 4> C /T C that is well above 
the critical first-order boundary ~ 1. 



IX. CONCLUSIONS 

We have calculated the strength of the electroweak 
phase transition in a supersymmetric model with four 
chiral generations. We find there is an intriguing region 
of parameter space, with fourth generation quarks heav- 
ier than about 300 GeV and the squark to quark mass 
ratio 1 < rriqi /m q i < 1.1, where <f> c /T c > 1. Within this 
region of parameter space we showed the Higgs can be 
easily heavier than the LEP bound of 115 GeV. 

This suggests that a viable model of electroweak baryo- 
genesis could indeed be a low energy supersymmetric 
model with a fourth generation of chiral fermions. What 
we have shown is the the strength of the first-order phase 
transition can be large enough to prevent the washout of 
a baryon asymmetry. This model also has several new 
sources of CP violation, ubiquitous in low energy su- 



persymmetry, that could be used to satisfy Sakharov's 
CP violation criteria. Examples of sufficient CP viola- 
tion that have been employed in other supersymmetric 
electroweak baryogenesis scenarios [5J DZ3 E] include the 
phase of the Higgsino mass parameter /i as well as the 
gaugino mass parameters M\^. 

It is coincidental that the region of parameter 
space where the first-order transition is strong enough 
combined with obtaining a large enough Higgs mass 
(taking tan/3 = 1) happens to be just beyond the 
current Tevatron direct search bounds 023 EH] . If we are 
lucky, the Tevatron could begin to see evidence for new 
physics in the form of both an extra chiral generation 
as well as superpartners in the very near future. The 
LHC, however, can easily cover this parameter space. 
Indeed, the mechanism to enhance the first-order phase 
transition described here is expected to be found or 
ruled out with only modest amount of data from the 
LHC. 



Note added: As this paper was being completed, 
Ref. [33] appeared, speculating that the electroweak 
phase transition could be enhanced with a fourth gener- 
ation without supersymmetry. Unfortunately, this does 
not work, as we show in Fig. [3] where just adding fermions 
actually decreases 4> C /T C because of the effects of the 
fermions on the zero-temperature effective potential. 



7 




(a) 



350 400 450 

rry [GeV] 




(b) 



350 400 450 

77v [GeV] 



FIG. 7: Same as Fig. [6] but we take m t i/m b i — 1.2, as favored by electroweak precision data. The basic shape and size of the 
region is the same, illustrating that our results are not particularly sensitive to the heavy fourth generation quark mass ratios. 



Acknowledgments hospitality where part of this work was done. This work 

was supported in part by the Department of Energy un- 
We thank M. Ibe, D. Reeb, T. Roy, G. Servant, T. Tait, der contract DE-FG02-96ER40969. 
and C. Wagner for discussions. GDK thanks CERN for 



[1] A. D. Sakharov, Pisma Zh. Eksp. Teor. Fiz. 5 
(1967) 32 [JETP Lett. 5 (1967 SOPUA,34,392-393.1991 
UFNAA,161,61-64.1991) 24]. 

[2] C. Jarlskog, Phys. Rev. Lett. 55, 1039 (1985). 

[3] G. W. Anderson and L. J. Hall, Phys. Rev. D 45, 2685 

(1992) . 

[4] M. Dine, P. Huet and R. L. . Singleton, Nucl. Phys. B 

375, 625 (1992). 
[5] P. Arnold and O. Espinosa, Phys. Rev. D 47, 3546 

(1993) [Erra tum-ibid. D 50, 6662 (1994)] [arXiv:hep-| 
ph/9212235 . 

[6] A. G. Cohen, D. B. Kaplan an d A. E. Nelson, Ann. Re v. 

Nucl. Part. Sci. 43, 27 (1993) [arXiv:hep-ph/9302210| . 
[7] M. Quiros, arXiv:hep-ph/9901312 

[8] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, 

Phys. Lett. B 155, 36 (1985). 
[9] G. F. Giudice, Phys. Rev. D 45, 3177 (1992). 
[10] J. R. Espinosa, M. Quiros and F. Zwirner, Phys. Lett. 

B 307, 106 (1993) |arXiv:hep-ph/9303317| ; A. Brignole, 

J. R. Espinosa, M. Quiros and F. Zwirner, Phys. Lett. B 

324, 181 (1994) [arXiv:hep-ph79 312296 . 
[11] M. S. Carena, M. Quir os and C. E. M. Wagn er, Phys. 

Lett. B 380, 81 (1996) | arXiv:hep-ph/9603420| . 
[12] D. E. Morrissey, "LHC Probes of (MSSM) Baryogenesis" , 

talk given at Physics of the LHC, KITP. 
[13] M. Pietroni, Nucl. Phys. B 402, 27 (1993) arXivihep- 



ph/9207227] , 

[14] J. Kang, P. Langacker, T. j. Li and T. Liu, Phys. Rev. 

Lett. 94, 061801 (2005) [arXiv:hep-ph/0402086l . 
[15] A. Menon, D. E. Morrisse y and C. E. M. Wagner , Phys. 

Rev. D 70, 035005 (2004) [arXiv:hep-ph/0404184] . 
[16] S. W. Ham, S. K. OH, C. M. Ki m, E. J. Yoo and D. Son , 

Phys. Rev. D 70, 075001 (2004) [arXiv:hep-ph/ 0406062] . 
[17] M. S. Carena, A. Megevand, M. Quiros and 

C. E. M. W agner, Nucl. Phys. B 716, 319 (2005) 

[arXiv:hep-ph/0410352] . 
[18] K. Funakubo, S. Tao and F. Toyoda, Pr og. Theor. Phys. 

114, 369 (2005) [arXiv:hep-ph/0501052| . 
[19] S. J. Huber, T. Konstandin, T. Prokopec and 

M. G. Schmidt, Nuc l. Phys. B 757, 172 (2006) 

|arXiv:hep-ph/0 606298 ; S. J. Huber, T. Konstandin, 

T. Prokope c and M. G. Schmidt, Nucl. Phys. A 785, 

206 (2007) [arXiv:hep- ph/06080i7l . 
[20] J. Shu, T. M. P. Tait and C. E. M. Wagner Phys. Rev. 

D 75, 063510 (2007) [arXiv:hep-ph/0610375] . 
[21] S. Profumo, M. J. Ramsey-Musolf and G. Shaughnessy, 

JHEP 0708, 010 (2007) [arXiv:0705.2425| [hep-ph]]. 
[22] C. Grojean, G. Se rvant and J. D. Wells Phys. Rev. D 

71, 036001 (2005) [arXiv:hep-ph/0407019| . 
[23] S. W. Ham, S. K. Oh and D. Son, Phys. Rev. D 71, 

015001 (2005) [arXiv:hep-ph/0411012j . 
[24] G. D. Kribs, T. Plehn, M. Spannowsky and T. M. P. Tait, 



8 



Phys. Rev. D 76, 075016 (2007) [arXl^ :0706.3718 [hep- 
ph]]. 

[25] M. Fukugita and T. Yanagida, Phys. Rev. D 42, 1285 
(1990); A. E. Nelson and S. M. Barr, Phys. Lett. B 246, 
141 (1990). 

[26] We thank M. Ibe for discussions on this point. 
[27] J. F. Gunion, D. W. McKay and H. Pois , Phys. Lett. B 
334, 339 (1994) [ arXiv:hep -ph/9406249 ; J. F. Gunion, 
D. W. McKay and H. Po is, Phys. Rev. D 53, 1616 (1996) 
[arXiv:hep-ph /9507323| ; M. S. Carena, H. E. Haber 
and C. E. M. Wagner Nucl. Phys. B 472, 55 (1996) 
|arXiv:hep-ph/9512446] . 
[28] D. M. Pierce, J. A. Bagger, K. T. Ma tchev and 
R. j. Zhang Nucl. Phys. B 491, 3 (1997) | arXiv:hep- | 
ph/9606211 



3] T. Aaltonen et al. [CDF Collaboration], |arXiv:0801. 3877 
[hep-ex] . 

[30] DO collaboration, http : //www-dO . f nal . gov/Run2Physics/ 
WWW/results/prelim/HIGGS/H21/H21F0_supplement . eps 

[31] S. P. Martin, |arXiv:hep- ph/970935"6l 

[32] M. S. Carena, H. E. Haber, S. Heinemeyer, W. Hollik, 

C. E. M. Wagne r and G. Weiglei n, Nucl. Phys. B 580, 

29 (2000) |arXiv:hep-ph/0001002| . 
[33] M. S. Carena, M. Quiros an dC. E. M. Wa gner, Nucl. 

Phys. B 524, 3 (1998) [arXiv:hep-ph/9710401| . 
[34] M. Laine and K. Rummukainen, Nucl. Phys. Proc. Suppl. 

73, 180 (1999) arXiv:hep-lat/9809045 . 
[35] We thank T. Roy for nudging us to consider counting 

differently. 

[36] W. S. Hou, |arXiv:0803.1234| [hep-ph]. 



